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We show that and how ultra-cold atoms in an accelerated two-band lattice are a controlled real- 
ization of Landau-Zener-Stiickelberg interferometry. 



I. INTRODUCTION 

Two-level systems subject to a strong periodic driv- 
ing appear in a large variety of quantum mechanical sys- 
tems and their study has a long history. They are nat- 
urally used in atomic physics when an atom is coupled 
to a laser-field pQ, in solid state physics when describ- 
ing superconducting qubits [2] or as effective models [3]. 
One consequence of the strong driving is the possibility 
for the two-level system to undergo a sequence of transi- 
tions. Each transition can be seen as an effective beam 
splitter and the coherent passage through several transi- 
tions leads to an accumulation of phases and interference 
effects known as Stiickelberg oscillations (see [1] and ref- 
erences therein). 

Recently, Stiickelberg oscillations have been observed 
experimentally for ultra-cold atoms in accelerated optical 
lattices [SJjBJ- This opens the route for very detailed stud- 
ies of Stiickelberg interferometry with cold gases. The 
high degree of control in these systems [3 |8] allows to 
explore the strong sensitivity of the phase between inter- 
band transitions on the band structure. The main goal of 
the present paper is to establish explicitly the connection 
between Stiickelberg interferometry and interband tran- 
sitions in optical lattices. Additionally, we provide simple 
analytical formulae for the interband dynamics and com- 
pute interference patterns, i.e. contour plots of transition 
probabilities, for realizations with a single optical lattice 
as experimentally used in [5] and superlattices as in [5J. 
The outline is as follows. In section [IT] we are going 
to obtain the Landau-Zener-Stiickelberg (LZS) Hamil- 
tonian from a two-band model for ultra-cold atoms in 
accelerated optical lattices and study the dynamics of 
the interband transitions using a systematic expansion 
in section |III| We will then compute the transition prob- 
abilities using degenerate perturbation theory and com- 
pare the results to numerical simulations. This will be 
followed by predictions for interference patterns in real- 
izations with a single optical lattice and superlattices. 
We will close with a short summary. 



II. COLD ATOM REALIZATION OF THE 
LANDAU-ZENER-STUCKELBERG 
HAMILTONIAN 

A quantum mechanical two-level system with energy 
bias £o under strong periodic driving with amplitude 
A and frequency uj is modeled by the Landau-Zener- 



Stiickelberg Hamiltonian [4] 

^lzs - -2 ^ At -eo-Atmutj- (1) 

Here, A^ denotes the tunneling amplitude between the 
two levels. In the present section we are going to show ex- 
plicitly how this Hamiltonian can be realized with ultra- 
cold atoms in accelerated optical lattices. 

Using ultra-cold atoms in optical lattices it is possi- 
ble to create tilted non-interacting two-band systems [5J . 
We make all parameters of the Hamiltonian dimension- 
less by measuring them in units of recoil energies E Tec = 
h 2 k\/ (2m), where kr, is the wave vector of the laser cre- 
ating the optical lattice and m the mass of the atoms 
(we set h = 1). The dimensionless force F is obtained by 
multiplication of the physical force with the lattice con- 
stant and dividing by E Tec . The appropriate dimension- 
less two-band Hamiltonian obtained from an expansion 
in Wannier functions of the full problem [9l [10] reads 

n 2B = [( lF - fhUi - y (oJ+iOj + h - c -) ( 2 ) 

+ (IF + f - y + h.c.) + FC (bj ai + h.c.) . 

The operator a/ (at) annihilates (creates) a particle at 
site I and bi (b t ) in the upper band. The bands are sep- 
arated by a bandgap A and the whole lattice is tilted 
by on-site energies IF. The hopping amplitudes between 
neighbouring sites in band a, b are denoted by J a > 
and Jb < 0. The single-particle coupling of the bands 
is proportional to the external Stark force F via CqF 
with a coupling constant Co depending on the optical 
lattice but usually of order C -0.2 [TD]. We 
take the external force F as a free parameter. The pa- 
rameters of the Hamiltonian are directly computed from 
Wannier functions wf' b (x), which are maximally local- 
ized states centered around the l-th lattice well. They 
can be computed for realizations with a single opti- 
cal lattice [TU] V(x) = V cos(a;) or with a superlattice 
V(x) — Vi cos(a;) + V% cos(2x + 4>) with a possible addi- 
tional phase cj) between the lattices [TT] [T2] . The param- 
eters are then given by 

J a , 6 = J w a ' h {x)V(x)w a ' b {x)dx (3a) 

C = / w a (x) -x-w b {x)dx. (3b) 
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Here, w?' (x) denote the Wannier functions for the low- 
est (a) and first excited band (b) that are computed from 
the Bloch functions to the potential V{x). The band gap 
A is the difference between the average energy of these 
two lowest Bloch bands. Inter-particle interactions can 
be neglected under specific experimental conditions [5] 
and the case of weak interactions can be treated pertur- 
batively, extending the discussion in [T3J [T3] . 

To obtain a driven time-dependent two-level system 
similar to eq. 0, we change to the interaction-picture 
with respect to the external force [T5]. This removes 



tilted optical lattice Landau-Zener-Stiickelberg 



the tilt J2i lF a } a i an d replaces 



i+i u i 



a, — > 



iFtJ 



then 



(and likewise for 6[ +1 6 ; ). The Hamiltonian is 
time-dependent with a periodicity of Tb = 2tt/F. In- 
troducing Fourier components a(k) = J^i e- llk ai and 
b(k) — X), e" fy, we obtain the following periodic two- 
level Hamiltonian [TBI EZ] 



-f - J a cos(fc + Ft) 
C F 



A 
2 



C F 
J b cos(k + Ft) 



(4) 



The finite distance between the two levels, a sinusoidal 
driving and a constant coupling as in the LZS Hamilto- 
nian, eq. 0, are already present now. To make the con- 
nection completely transparent, we add a periodic shift 
of the energy zero. The final result is then 



H = — 



1 / A + 2 J cos(fc + Ft) -2C Q F 
2\ -2C F -A - 2Jcos(fc + Ft) 

1 
1 



+ ( Ja + Jb) COs(k + Ft) 



(5) 



where J = J a — J&. Shifting the time zero as t — > t—k/F— 
7r/2F, we arrive exactly at the form of eq. ([1]). The first 
part of the Hamiltonian eq. ^ is the ultra-cold atom re- 
alization of the Landau-Zener-Stiickelberg Hamiltonian 
and the second part reflects a time-dependent shift of the 
zero energy point, which does not concern interference 
effects between different phases. To simulate the LZS 
Hamiltonian and to perform interferometry one can use 
the control offered by ultra-cold atom systems in acceler- 
ated optical lattices. The role of the energy bias is taken 
by the average band gap between the two Bloch bands, 
the driving amplitude is realized as the difference in hop- 
ping strengths and the driving frequency is the Bloch 
frequency lob = F in our units. Finally, the tunneling 
amplitude in the ultra-cold atom realization is propor- 
tional to the external Stark Force. This correspondence 
between the two realizations is summarized in table |U 
Please note that not all parameters in the cold atom re- 
alization can be varied independently since the driving 
frequency ljb — F and the interband coupling CqF both 
depend on the external force. Figure [T] presents a typical 
LZS interferometric pattern (see [3] for many examples 
and a summary of different experimental results) . Shown 
are the transition probabilities, i.e. the long-time aver- 
age (t ^> Tb) of the occupation of the upper band, when 
varying the 'driving strength' J and the 'level splitting' 



band gap A 
Bloch frcq. ojb — F 
band coupling 2CoF 
hopping J = J a — J b 



level splitting eo 
driving freq. ui 
level coupling — At/2 
field strength A 



TABLE I: Analogy between tilted optical lattices and the 
Landau-Zener-Stiickelberg Hamiltonian. Note that the band 
coupling and the driving (or Bloch-) frequency are not inde- 
pendent in tilted optical lattices. 
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band gap A 

FIG. 1: (Color online) Landau-Zener-Stiickelberg interferom- 
etry. Shown is the long-time average of the occupation of the 
upper band as a function of the band gap A (corresponding to 
the level separation in atomic systems) and the hopping dif- 
ference between the bands J = J a — Jb (corresponding to the 
driving amplitude). Parameters: Co = —0.15 and F = 1.0. 



A according to eq. (21) to be derived below. To see how 



these interference patterns arise in systems of cold atoms, 
we will derive analytical expressions for the interband dy- 
namics. 



III. STUCKELBERG INTERFEROMETRY 
WITH ULTRA-COLD ATOMS 



A. Dynamics of the interband transitions 

To study the dynamics of the interband transition we 
have to solve the time-dependent Schrodinger equation. 
However, it can be shown to be equivalent to the Hill 
equation and exact analytical solutions in closed form are 
not possible. To obtain approximate solutions in a sys- 
tematic fashion, we will use the Magnus expansion [18j . 
It is useful apply the following transformation in order to 
obtain a purely off-diagonal Schrodinger equation [16U17] 
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a(k, t) = a(k, t) e - <A - t / 2 - <J » K «-( fc +«')d*' (6a) 
6(k, t) = b{k, t )e +lA - t / 2 - lJb /o co S (fc+Ft')dt' _ ( gb ) 

This removes the diagonal terms and we obtain an equiv- 
alent Schrodinger equation with a transformed Hamilto- 
nian Wit) for the transformed amplitudes 
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'dt 



a(k, t) 
b(k,t) 



/ CoFe-^A (a{k,ty 
[CoFe***,*) J \b(k,t) y 

(7) 

where (p(k, t) = A • t - (J/F)[sm(k + Ft) - sin(fc)] is the 
phase between the two Bloch bands and J = (J a — Jb). 
Where the pure existence of two energy bands allows phe- 
nomena like Rabi oscillations, it is the non-trivial phase 
difference for the time-evolution in both bands, related 
to J ^ 0, that gives rise to Stuckelberg oscillations and 
the complex interference phenomena. That is, the differ- 
ence in the curvature J = J a — Jb ^ allows to collect 
different phases during the time evolution in the upper 
or lower band and excludes exact analytical solutions of 
the Schrodinger equation in closed form [TJ [H [T7] . 

The idea of the Magnus expansion is to express the 
time evolution operator as the exponential of an infinite 
series U(t) — exp \y~) 1 Q n (t)] , with each term contain- 
ing an increasing number (n — 1) of nested different time 
commutators [18]. All orders of the Magnus expansion 
could be summed in a recent work for a new derivation 
of the transition probability in the Landau-Zener prob- 
lem [19] . The first two terms of this series read explicitly 
(valid in both Schrodinger and interaction picture |20j ) 



J o 

I rt rti 



(8a) 



*M*) = o / d<1 / d<2 Mt^Hih)} (8b) 
* Jo Jo 

To compute the terms in the Magnus expansion for our 
time-dependent problem, eq. ([7]), explicitly, we have to 
integrate the Hamiltonian and commutators of it over 
time. Without loss of generality, we restrict our discus- 
sion to k = 0, which can always be achieved by shifting 
the time zero. We will need the integral 



x(t) 



>df 



,i(A-t'-(J/F)sinFt') d1 J_ /gN 



Using the the generating function of the Bessel function 
exp[w(r — 1/t)] = J2nez Jn(2u)r n for r = e~ %Ft , we can 
write the sin in the exponent as a sum over Bessel func- 

tions X (t) = EnezMJ/F) /*df e i < A -» F >*'. After inte- 
gration and minor manipulations, we obtain a closed ex- 
pression with the explicit time dependence for the highly 
oscillatory function 



;(t) = 2^J„(J/F) e 



sin (w n t/2) 



(10) 



with ui n = A — nF. With the explicit expression for the 
integral defining x(£), the first term in the Magnus expan- 
sion is given by f2i(t) = —iCqF^ ® X Q ) and the time evo- 
lution operator in first order reads correspondingly [11) 



Ui(t) 



cos(C F|x|) 
-i^sXsm(fi< F\x\ 



3 iar s^sin(C F|x| 
cos(C F| X |) 



This is the result in first order and the occupation of the 
upper band Pb(t) = J2k\b(k,t)\ 2 is given by 



Pb(t) 



2C Q F 



E 



MJ/F)e i 



,t/2 sin (u> n t/2) 



(12) 

This result captures resonant and non-resonant contri- 
butions of the interferometry in a single and explicit for- 
mula. Eq. ( 12 ) can be understood by treating the denom- 
<^ A = mF contained in x(t) sep- 
>0 = t. This condition A w mF 



mator zeros w m = 
arately using linx, 



corresponds to a resonant interband coupling and x(t) 
can be decomposed as 

X {t) = J m {J/F) t + 2^ UJ/F) e^' 2 SinKV2) . 

For large times the first term will be dominating and the 
overall short-time averaged occupation of the upper band 
shows large sinusoidal oscillations with unit amplitude 
Af b ICS (t) = sin 2 [V, J m { J /F)t]. The other high-frequency 
and non-resonant terms lead to small amplitude oscilla- 
tions on top of this overall resonant interband oscilla- 
tions. 

For the second order contribution we need the commu- 
tator (where a z is the diagonal Pauli matrix) 

[H(ti),n(t 2 )} = 2iC 2 F 2 a z su#(t a ) ~ <K*i)] (13) 

and integrate over it in time. One obtains i7 2 (*) = 
iC§F 2 o z T\)(t) where the required integral ip(t) = 

J^dti / ' 1 dt2 sin[0(t2) — 4>(ti)] can be computed by ap- 
plying the same expansion as above 



LOn l- 



n£2 



' sin 2 u n 



m£2 



t/2 



'[(m + n)Ft/2] 
(m + n)F 



)}■ 



(14) 



The time evolution operator can again be given ex- 
actly |11) and one finds for the occupation of the upper 
band 



Pb(t) = 



\ x *x + C 2 F^\ 



sin 2 2C F 



X*x + C 2 F 2 r 



(15) 

where we suppressed the time-dependence of the func- 
tions x(t) an d V'W- The eqs. (12) and (15) are the cen- 
tral results of this section. They provide good approxi- 
mations to the full interband dynamics in an explicit ex- 
pression. They furthermore contain the resonant as well 
as non-resonant contributions to the interband dynamics. 



4 



B. Average occupation of bands 

In order to determine the average occupation of the 
bands it is more useful to go back to the original two- 
band Hamiltonian eq. ^ instead of calculating the 
long-time average over the expressions eq. ( [l2| and 
eq. (15 1. Introducing the Wannier-Stark states [21] 

«» = Ll£7^-n(^/ f ) Q ! and Pn = J2leZ J l-n( J b/ F )bl 

into eq. pF, one obtains [T3] 



2B 



IF 



5^ C FJ m (J/F)(afa- m + h.c. 



meZ 



(16) 



where J — J a — Jb as above. This expression contains 
all relevant processes coupling the two bands as direct 
couplings weighted by Bessel functions J m (J/F). For 
J/F < 1 the dominant contribution is the on-site cou- 
pling between the bands proportional to Jq(J/F). Keep- 
ing only this dominant contribution, the Hamiltonian 
eq. ( 16 1 is a sum of independent two-level systems and 
can easily be diagonalized. The resulting occupation of 
the upper band contains the dominant part of the non- 



resonant oscillations already contained in eq. ( 12 ) and 
reads 



P (t) W » 



Sill 



A 2 



t/2 



(17) 



where Vq = C^FJ (J/F). This means that the non- 
resonant contribution to the averaged occupation of the 
upper band is given by 



Ph = 



1 



1/2 

r ^ ■ 

12C J (J/F). 



(18) 



However, in addition to this onsite coupling, the 
Hamiltonian eq. (16) allows a direct coupling of more 



remote sites. This becomes particularly important when- 
ever a site from the lower band and a site from the upper 
band are energetically degenerate. The corresponding 
resonance condition for two levels being separated by m 
sites is Aw mF, m £ N. We therefore apply degenerate 
perturbation theory 22^ to the Hamiltonian eq. ( 16 ) and 
obtain for a resonance of order m in second order the 
following effective two-level system [TTJ [53] 



\Vi- m ^\ 2 



+ 

i^l c 2-m 
V-rn 



V—r 



\Vi- 



i^l — m * 



(19) 



with ef = IF ± A/2 and Vi = C FJi(J/F). This and 
higher orders allow the computation of various observ- 
ables with high degree of precision [53] . For example, the 
resonance condition A = mF experiences a slight Stark 



shift and the corresponding condition in second order is 
given by 



A = mF — 2C 2 F 2 ^ 



A-iF ' 



(20) 



Unlike the usual LZS problem, the coupling between the 
bands CqF and the driving frequency ujb — F are not 
independent for atoms in optical lattices. This makes 
eq. ( [20] ) nonlinear and difficult to solve. However, it 
can be solved either numerically or by iteration. The 
uncorrected resonance position for a single optical lat- 
tice with Vo = 4 and a resonance of order 2 is given 
by F 2 = A/2 = 2.195. Solving eq. p0| numerically 



gives -Fjf T = 2.22067 which has a relative error of order 
10~ 5 when compared to the maximum of the resonance 
F = 2.22070 from numerical simulations of the full prob- 
lem eq. Q. In the same way, very high precision can be 
achieved by higher orders perturbation theory. 

The resulting resonant contribution to the average oc- 
cupation of the upper band has a Lorentzian shape [41,117] 
and the total transition probability is given by the non- 
resonant interband coupling and the different resonant 
contributions 



Pi 



1 



2 A 2 



4(V m /FA) 2 



(1/F-1/F m ) 2 +4(VWFA) 2 ' 

(21) 

The first term describes the direct force- induced coupling 
between the bands and is usually much smaller than 
unity. However, for F 3> A its contribution becomes 
important and grows as F~ 2 . The second part are the 
resonant contributions from different orders of resonance. 

to numerical simula- 



Let us compare this result eq. (21 



We change to the 
le external force and 



tions of the full Hamiltonian eq. 
interaction picture with respect to t 
impose periodic boundary conditions [TUIUS] ■ We take an 
initial state \tpo) which occupies only the lower band, and 
evolve it in time according to the Schrodinger equation 
idt\ip(t)) = H(t)\ip(t)). We compute the total occupa- 
tion of the upper band Pb(t) = (*0(t)|&J&/|^>(t)) and take 
the long time average. The result for a single optical lat- 
tice realization with Vq = 4 is shown in fig. [2] together 
with our result eq. (21). We observe very good agree- 



ment even on the logarithmic scale shown in the figure. 
The non-resonant interband coupling and the resonant 
contribution as well as the resonance positions are accu- 
rately reproduced. Only the asymmetry of the resonance 
peaks (in particular of the high order m = 4 resonance 
at l/F w 0.9) is not captured by our present analysis 
since the effective model of eq. ( fl9[ ) should be extended 
to describe well such more complex peak profiles. 



C. Interferometry with optical lattices 

We have established the general possibility to use ultra- 
cold atoms for Stiickelberg interferometry and have given 
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FIG. 2: (Color online) Longtime average of the occupation of 
the upper band from numerical simulations of the full problem 
(blue circles) and the theoretical prediction eq. (21 1 (dashed 
line). The insets show the resonances m = 2,3,4 on a linear 
scale. Parameters for a single optical lattice of depth Vb = 4: 
Co = -0.14, A = 4.39, and J = -0.682. 



analytical results for the probability of interband tran- 
sitions in the previous paragraphs. As already men- 
tioned, the important parameters hopping J and band 
gap A to build the Landau-Zener-Stiickelberg Hamilto- 
nian, eq. (JlJ, can, however, not be controlled indepen- 
dently in optical lattice systems. We therefore show in- 
terference patterns similar to figure [I] but with the ex- 
perimentally accessible parameters varied. That is, we 
vary the depth of the optical (super-)lattice (and possi- 
ble the phase <fi) and compute the Wannier functions for 
each value of the lattice parameters. From these we ob- 
tain the relevant parameters J, A, and Co according to 
eq. (J3j) . With them we obtain the transition probabilities 
at different forces from eq. (21 1. The results for realiza- 



tions with a single optical lattice and a superlattice are 
shown in fig. [3j We clearly observe resonances of differ- 
ent orders as the external force is varied. The resonance 
position changes nonlinearly with the lattice depth since 
the band gap is generally a not strictly linear function of 
the lattice depth (in both cases of a single lattice and a 
superlattice). Additionally, the importance of the non- 
resonant background becomes more important for larger 
forces, i.e. for small 1/F. 

As mentioned already, the three system parameters J, 
A, and Co all depend on the lattice depth Vq (or V2/V1) 
and cannot be varied independently. However, the situa- 



tion is slightly advantageous for superlattices, since there 
are two experimental parameters, the ratio of the lattice 
depths V2/V1 and the relative phase between the lattices 
</>, that can be altered. We therefore computed the Wan- 
nier functions and the system parameters for many dif- 
ferent combinations of these two parameters and show 
the resulting transition probability as a contour plot in 
fig. [4] for fixed external force F = 3. We observe a clear 
and broad resonance as a result of Stiickelberg interfer- 
ence when different parts of the wave function evolve in 
the different bands. The complicated shape of the res- 
onance is again a result of the nonlinear dependence of 
the system parameters on V2/V1 and <j>. Fig. [4] is an ex- 
plicit prediction for transition probabilities from multiple 
phase interference that should be observable with current 
experimental methods as in [5J. 



IV. SUMMARY 

We have made the connection between Stiickelberg 
interferometry and recent experiments with ultra-cold 
atoms. We showed explicitly how to obtain the Landau- 
Zener-Stiickelberg Hamiltonian with cold atoms in ac- 
celerated optical lattices. More specifically, we applied 
the Magnus expansion to obtain analytical expressions 
capturing various aspects of the complicated interband 
dynamics. The transition probabilities for different ex- 
perimental realizations with atomic quantum gases have 
been computed and should be experimentally accessible. 
We thus hope to have clarified some of the background of 
ongoing experiments and to stimulate further research us- 
ing the high control in state-of-the-art implementations. 
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